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Abstract 

This paper compares two methods for predicting 
transonic rotor noise for helicopters in hover and for- 
ward flight. Both methods rely on a computational 
fluid dynamics (CFD) solution as input to predict the 
acoustic near and far fields. For this work, the same 
full-potential rotor code has been used to compute the 
CFD solution for both acoustic methods. The first 
method employs the acoustic analogy as embodied in 
the Ffowcs Williams-Hawkings (FW-H ) equation, in- 
cluding the quadrupole term. The second method uses 
a rotating Kirchhoff formulation. Computed results 
from both methods are compared with one other and 
with experimental data for both hover and advancing 
rotor cases. The results are quite good for all cases 
tested. The sensitivity of both methods to CFD grid 
resolution and to the choice of the integration sur- 
face/volume is investigated. The computational re- 
quirements of both methods are comparable; in both 
cases these requirements are much less than the re- 
quirements for the CFD solution. 

Introduction 

Modern helicopter rotor designs require both high 
aerodynamic performance and, to an unprecedented 
degree, low noise radiation. Low noise radiation is par- 
ticularly important for civilian helicopters, which can 
only offer significant operational benefits if they can 
operate in and nearby heavily populated areas. Great 
progress has been achieved in the accurate and reliable 
prediction of helicopter rotor noise. A recent review 
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of helicopter noise prediction is given by Brentner and 
Farassat [1], 

Among the various helicopter rotor noise sources, 
impulsive noise is particularly annoying as well as dif- 
ficult to predict accurately. The first type of impul- 
sive noise, high-speed impulsive (HSI) noise, is caused 
by compressibility effects associated with high-speed 
blade motion. This HSI noise is generally associated 
with advancing tip Mach numbers above Mat = 0.85 
and with the appearance of shock waves in the flow 
Held around the rotor. A second type of impulsive 
noise comes from the interaction of the rotor blades 
with their vortical wake systems. Blade- vortex interac- 
tion (BVI) is a significant source of noise near civilian 
population centers because it commonly occurs when 
helicopters are descending for landings. During de- 
scent, the advancing tip Mach numbers are usually 
subcritical. In this paper, we concentrate on HSI noise; 
however, our methods are readily applicable to BVI as 
well. We did not consider BVI noise so that we could 
avoid the errors associated with the wake prediction 
and, thus, simplify the comparison. 

The aerodynamics of HSI noise can be computed 
accurately with nonlinear computational fluid dynam- 
ics (CFD) methods based on full-potential or Euler 
and Navier-Stokes equations (e.g., [2] and [3]). How- 
ever, CFD codes are not practical for calculations far 
(i.e. , beyond 2 to 3 rotor radii) from the rotor blades. 
Thus, near-field CFD solutions have been coupled with 
more efficient computational methods to propagate the 
acoustic signals to the far field. 

The Lighthill acoustic acoustic analogy [4, 5], as ex- 
pressed by the Ffowcs Williams-Hawkings (FW-H ) 
equation [6], has been used for nearly 20 years to pre- 
dict noise from helicopter rotors. The acoustic anal- 
ogy predicts discrete-frequency noise quite well, but 
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cannot predict HSI noise if the quadrupole term is ig- 
nored. Thus, acoustic analogy prediction of HSI noise 
was forced to wait until CFD methods could predict 
the transonic flow held around rotor blades moving at 
high speed. 

After advancements in CFD capability enabled ac- 
curate computation of the three-dimensional unsteady 
how held, another prediction methodology, the Kirch- 
hoff method, became viable as well. The Kirchhoff 
method assumes that all sound sources are computed 
correctly by a CFD code and that the acoustic waves 
are propagated to a surface that surrounds all the 
sources. When the pressure, the temporal derivative 
of pressure, and the normal derivative of pressure on 
the surface are known, the Kirchhoff formulation can 
determine the acoustic pressure outside the surface via 
linear wave propagation. 

The main purpose of this paper is to compare and 
contrast these two acoustic prediction methodologies. 
This comparison is useful because both methods re- 
quire essentially the same CFD calculation to provide 
the input data; however, whether one method has an 
advantage in terms of efficiency, accuracy, and robust- 
ness over the other is still unclear. The CFD code 
used here, as well as both acoustic prediction methods, 
is implemented in the TiltRotor Aeroacoustic Code 
(TRAC) [7] at NASA Langley Research Center; we 
expect that the results of this paper will be useful to 
future users of this code system. 


Acoustic Prediction Methods 

In this section, a brief description of the two acous- 
tic prediction methodologies is given, along with some 
details of the codes used in this study. The acous- 
tic codes that are compared in this work have similar 
numerical algorithms and require approximately the 
same computer resources. The CFD code used in this 
work will also be discussed. 


Acoustic Analogy 

The FW-H equation [6] is the most general form of 
the Lighthill acoustic analogy and is appropriate for 
predicting the noise generated by the complex mo- 
tion of helicopter rotors. In differential form, the 
FW-H equation is given by the following inhomoge- 
neous wave equation: 


n 2 p'(x,t) 




d 2 

dx;dx , 


■[TijH(f)] 


(1) 


where p'(x,t) is the acoustic pressure and the three 
source terms on the right-hand side are known as the 


thickness, loading, and quadrupole source terms, re- 
spectively. The rotor blade is defined by the equation 
/ = 0, v n is the local normal velocity of the surface, 
li are the components of the local force vector exerted 
by the surface on the fluid, and Tij is the Lighthill 
stress tensor. The FW-H equation is valid in the en- 
tire unbounded space; hence, a formal solution may 
be obtained by using the free-space Green’s function 
8(g) /4irr. 

Note in equation (1) that the thickness and load- 
ing source terms are surface distributions of sources 
(indicated by the presence of the Dirac delta function 
8(f)). These source terms have been used for sev- 
eral years in rotor noise prediction because they ac- 
count for most of the acoustic signal when the flow 
Held is not transonic and they do not require knowl- 
edge of the field around the blade (although the ac- 
curate determination of the blade-surface pressure is 
still challenging). The quadrupole source, on the other 
hand, is a volume distribution of sources (indicated by 
the Heaviside function H(f)). The importance of the 
quadrupole term has long been recognized [8, 9]. How- 
ever, the quadrupole source has often been neglected 
in rotor noise prediction because of the computational 
demands of computing the flow field with sufficient 
accuracy and integrating over a volume in the acous- 
tic prediction. Recently, Brentner and Holland [10] 
developed an efficient and robust method to predict 
HSI noise for a hovering rotor by using a far-field ap- 
proximation of the quadrupole source. Brentner [11] 
extended this method for advancing rotors in the code 
now called WOPWOP+. The WOPWOP+ code is 
utilized in this work for all acoustic analogy predic- 
tions. 

WOPWOP+ 

The original WOPWOP computer code, developed 
by Brentner [12], is used to solve the FW-H equation 
for near- and far-field noise. The WOPWOP code 
uses Farassat’s formulation 1A [13], which models 
only the thickness and loading noise terms in the 
FW-H equation. The WOPWOP code requires time- 
accurate rotor surface pressures and their temporal 
derivatives for the loading noise prediction. Realis- 
tic helicopter blade motions can also be prescribed as 
inputs. The WOPWOP code is widely used in the 
rotorcraft community, which includes all major U.S. 
rotorcraft companies. 

A new far-field approximation to the FW-H 
quadrupole source, which is exact for an in-plane far- 
field observer, has recently been added to the WOP- 
WOP code. (Since HSI noise radiates primarily in 
the rotor plane, this approximation is also expected 
to work well for observers outside the rotor plane.) 
The new version, called WOPWOP+ [10, 11], uti- 
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lizes preintegration of the quadrupole volume source 
in the direction normal to the rotor disk to reduce 
the computational time required for the quadrupole 
noise calculation. This far-held approximation leads 
to quadrupole integrals of the same form as thick- 
ness and loading noise (i.e. , retarded-time surface in- 
tegrals). Quadrupole noise prediction is carried out in 
two parts: a preprocessing stage in which the previ- 
ously computed how held is integrated in the direction 
normal to the rotor disk and a noise computation stage 
in which quadrupole surface integrals are evaluated for 
a particular observer position. A retarded-time formu- 
lation is used for the quadrupole noise computation; 
thus, subsonic quadrupole source panel motion is re- 
quired. The quadrupole noise computation is robust 
and requires computer resources comparable to those 
required for thickness and loading noise prediction. 


Kirchhoff Formulation 

The second prediction method for far-held heli- 
copter acoustics is based on Kirchhoff ’s theorem. Pres- 
sures and pressure derivatives from the near-held CFD 
computation are utilized on a surface that encloses 
the rotor blades. The derivation of the Kirchhoff for- 
mulation depends on the assumption that only linear 
acoustic propagation occurs outside the Kirchhoff inte- 
gration surface and all acoustic sources and nonlinear 
effects are contained within the surface. Hence, the 
sound speed is constant for the acoustic wave propa- 
gation from the Kirchhoff surface to the far held. The 
most important nonlinearities are captured, however, 
because the CFD solver computes the complete non- 
linear acoustic held inside the Kirchhoff surface. The 
method was hrst demonstrated in rotorcraft aeroa- 
coustics for two-dimensional cases and was later ap- 
plied in three-dimensional cases. A general review 
of the use of the Kirchhoff method in computational 
acoustics is given by Lyrintzis [14]. 

Farassat and Myers [15] have derived a general 
Kirchhoff formula that gives the pressure signal of an 
observer in the stationary coordinate system as a func- 
tion of the surface integral over surface S of the pres- 
sure, the normal derivative of the pressure, and the 
time derivative of the pressure, where surface S is 
moving at an arbitrary speed. The formulation was 
derived for a deformable surface; however we use the 
rigid-surface version as follows: 


p(X*,t*) 
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where 

Ei = {Ml - 1 )p n + M n M t ■ V 2 p - a^M n p (3) 

a- 1 

+ ^ ) [(»r - Mn ~ n M )p+ (cos 9 - M n )p] 

+ ^ y [Mr (cos 6 - M n )p] 

^ = { \ZMrf {ZOSe ~ Mn) (4) 

The dot over M and n denotes the source time r* 

derivative. Here, (X*,t*) and (X, r) represent the ob- 
server position and observer time, and the source posi- 
tion and source time, respectively. The subscript r* in 
equation (2) indicates that the integrand is evaluated 
at the source emission time r* , which is the solution 
of the equation 

T-U + |X* -X(r)|/ aoo =0 (5) 

In addition, the following definitions are introduced: 

M r = M • f , M n = M • n, 
h r = nr, um = n • M ( 6 ) 

The formulation in equation (2) is valid when the 
observer is stationary and the surface is moving with 
arbitrary velocity. The case in which the observer is 
moving along with the rotor (e.g., rotor in a wind tun- 
nel with flow) can also be modeled simply by moving 
the observer position for each point in the time his- 
tory. A limitation of the formula is that singularities 
appear for supersonically moving surfaces. Farassat 
and Myers [16] have derived a new formula to avoid 
this problem; however the new formula has not yet 
been applied in practical problems. 

Two approaches have been used to implement 
the Kirchhoff integral for three-dimensional rotorcraft 
noise prediction. They differ as to whether or not the 
Kirchhoff surface rotates with the helicopter blades. 
The Hrst approach (e.g., [17] and [18]) employs a ro- 
tating Kirchhoff surface. This approach has an ad- 
vantage in that the Kirchhoff integrals use the same 
computational mesh as the CFD calculations; there- 
fore, interpolations onto the Kirchhoff surface are not 
required. One disadvantage of the rotating Kirchhoff 
implementation, however, is that parts of the surface 
sufficiently far from the center of rotation will poten- 
tially have supersonic velocity. The requirement for 
subsonic source motion imposes a limitation on the 
radial extent of the Kirchhoff surface for rotors whose 
tip Mach numbers approach unity. 

The second implementation approach (e.g., [19] and 
[20]) employs a nonrotating Kirchhoff surface that 
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completely surrounds the spinning rotor blades. In- 
terpolation of the input data is required to transfer 
the rotating CFD solution to the nonrotating Kirch- 
hoff surface; however, because the surface only trans- 
lates with the rotor, sonic velocities are avoided. Re- 
cently, Strawn et al. [21] compared the rotating and 
nonrotating Kirchhoff formulations for both HSI and 
BVI noise. Both methods yielded similar results for 
the cases considered. Note that high-accuracy CFD 
information (usually generated by an Euler code) out 
to the nonrotating Kirchhoff surface is needed. This 
requirement could increase the CFD needs of the non- 
rotating method. In this paper, the rotating Kirch- 
hoff implementation is used because the full-potential 
solver computes the solution only relatively close to 
the blade. In addition, the algorithmic details of the 
rotating Kirchhoff method are quite similar to those 
of the the acoustic analogy code. 

RKIR 

The numerical implementation of the Kirchhoff for- 
mula of Farassat and Myers for a rotating surface, 
known as RKIR, is described in reference [17]. The 
RKIR code requires an external source of data for pres- 
sure and pressure derivatives on a surface. In addi- 
tion to the time derivatives of pressure, the Kirchhoff 
scheme also requires spatial gradients of pressure. The 
computation of the temporal and spatial derivatives of 
the pressure are performed as a subroutine in the CFD 
flow solver. The Kirchhoff surface is typically located 
1 to 2 chord lengths away from the rotor surface and 
surrounds the rotor blades. The surface has the ap- 
pearance of a cylinder with a tip-end cap. The root 
end of the Kirchhoff surface is not included because 
its contribution is negligible. The Kirchhoff surface 
must be chosen to have only subsonic source motion 
for the reasons described previously. This restriction 
is the same as in the WOPWOP+ code. The RKIR 
code was recently modified [22] to include the same 
blade motions as WOPWOP (i.e. , pitching, lead lag, 
and feathering) and to make the input and output of 
RKIR as close to that of WOPWOP as possible. This 
version, which is part of the NASA Langley Research 
Center TRAC system [7], was used in this study. 

CFD Model 

A version of the full-potential rotor code FPR [2] 
was used to provide the flow-held computation in this 
work because it provides a fast and accurate CFD solu- 
tion, is widely accepted by the U.S. rotorcraft indus- 
try, and has recently been coupled to both acoustic 
codes exercised in this work. The FPR code uses a 
finite-difference scheme to compute three-dimensional 
unsteady transonic flows around the rotor blades. This 
code is fast in comparison with solutions of the Euler 


and Navier-Stokes equations and still gives a realis- 
tic representation of the physical flow field. The full- 
potential assumptions of inviscid and irrotational flow 
are appropriate for a transonic rotor. The flow field is 
dominated by inviscid effects, and the shocks are not 
strong enough to generate significant vorticity. 

The code was extended to accept wake information 
to account properly for rotor inflow and blade-vortex 
interactions in the version known as FPRBVI [23]. 
This version was used in the present study, although 
the wake and BVI capabilities were not utilized. All 
flow-field calculations computed here are for nonlift- 
ing rotors. The FPRBVI code has an internal alge- 
braic grid generator that is different from that in the 
original FPR code. This internal grid generator was 
used for all cases presented here. The near-field solu- 
tions from FPRBVI should provide accurate informa- 
tion for far-field calculations with both acoustic meth- 
ods. The same FPRBVI solution was used as input to 
both WOPWOP+ and RKIR. 

Results 

Acoustic predictions were made for three different 
model rotors operating both in hover and forward- 
flight test conditions. Two different grids were used 
in the FPRBVI computations. The first grid was a 
standard O grid with 80 points in the wrap-around di- 
rection, 36 points in the spanwise direction (27 points 
on the blade), and 24 points in the direction normal 
to the blade surface. A second, finer grid was also 
used, which had 150 points in the wrap-around direc- 
tion, 50 points in the spanwise direction (34 points 
on the blade), and 50 points in the normal direction. 
The computations were performed by first running 
FPRBVI in a quasi-steady mode for 400 time steps and 
then continuing in an unsteady mode. The unsteady 
solution was superior to the quasi-steady calculation 
for the hover cases, so this solution was used for both 
the hover and forward-flight calculations. The acous- 
tic codes, however, only used the CFD solution from 
one time step for the hover predictions. For the stan- 
dard grid, an unsteady time step of 0.125° azimuth was 
used; for the fine grid, a smaller time step of 0.0625° 
azimuth was required. The standard-grid computation 
required approximately 0.5 hr on a Cray C90, whereas 
the fine grid required approximately 5 hr. 

A comparison of computation time required by the 
acoustic codes is somewhat tentative because both of 
these codes are relatively immature. Nevertheless, we 
can compare of the computational resources needed by 
the current version of each code. The present codes de- 
pend upon file input / output (I/O) for data storage and 
manipulation during the prediction. File I/O can be as 
much as 25% or more of the total computational time 
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required for the acoustic predictions. For the case of a 
hovering rotor with the standard CFD grid, both codes 
require nominally 70 CPU sec per observer to run on 
an HP 9000 Series 735-99 workstation. For the fine 
grid, the WOPWOP+ code requires approximately 8 
seconds more while the RKIR code needed an addi- 
tional 39 seconds. The RKIR code required more time 
because the control surface for the fine grid contained 
nearly twice as many points; however, the source grid 
generated by the WOPWOP+ preprocessor did not 
change in size. For forward-flight, computations with 
the standard grid, WOPWOP+ required 791 sec to 
run (691 sec for the preprocessor and approximately 
100 sec for the acoustic computation) while the RKIR 
code only needed 520 sec. Much of the time expended 
by the WOPWOP+ preprocessor was associated with 
file I/O. The preprocessor is only run once per operat- 
ing condition; thus, if calculations are made for more 
that one observer location the WOPWOP+ computa- 
tion is faster. The processing time required to compute 
the temporal and spatial derivatives of pressure in the 
FPRBVI code is not accounted for in the Kirchhoff 
times. It is important to note that the computational 
time necessary to perform the acoustic predictions was 
substantially less than that required to compute the 
flow field. 

For the Kirchhoff predictions, a choice of the Kirch- 
hoff surface location is required. Several surface loca- 
tions, which correspond to grid lines, were evaluated 
in this work. Two surfaces, identified as the “large 
cylinder” and the “small cylinder,” respectively, are 
presented in this paper. For the standard grid, the 
large cylinder corresponds to the 21st. grid line normal 
to the blade surface; the small cylinder corresponds to 
the 18tli grid line normal to the blade surface. The 
tip-end cap for both surfaces was located on the 31st. 
spanwise grid plane. The location of these surfaces, 
shown for a standard grid in Fig. 1, was nominally 
maintained for all grid and rotor-blade combinations. 
The volume integration in WOPWOP+ used approx- 
imately the same volume contained by the Kirchhoff 
surfaces in Fig. 1. 

Sensitivity Analysis 

To determine the sensitivity and robustness of each 
prediction method, comparisons were made between 
the size of the volume around the rotor blade (i.e., the 
control surface position in the case of the Kirchhoff 
code) and the resolution of the CFD grid. These com- 
parisons were made for most of the cases presented in 
this paper; however, for brevity we present only the 
results for one rotor and operating condition. 

The first, comparison was for a. hovering rotor oper- 
ating at. a. high tip speed. The results of an experiment 



Figure 1. Location of the Kirchhoff surfaces (standard 
UH-1H grid). 

conducted by Purcell [24] are used for comparison. 
The rotor is a. one-seventh scale UH-1H model main 
rotor with straight., untwisted blades and an NACA 
0012 airfoil section. The rotor radius is 1.045 m with 
a. chord of 7.62 cm; thus, the blade aspect, ratio is 
13.7. The test was conducted with very little lift pro- 
duced by the rotor. Comparisons with predicted and 
measured acoustic signals are made for an in-plane mi- 
crophone located 3.09 rotor radii from the rotor hub. 
The experiment, was repeated for several Mach num- 
bers. Three hover Mach numbers, Mu = 0.70, 0.88, 
and 0.90, are used in this paper to compare the pre- 
dicted and measured results. 

The hover case with Mjj = 0.88 was chosen for the 
sensitivity comparison because the flow field is super- 
critical but not delocalized. This case should test the 
two acoustic codes because transonic effects are strong; 
however, this case still falls within the range of validity 
of WOPWOP+ and RKIR because supersonic source 
motion is not expected to be required. In Figs. 2 and 
3, a. comparison between the predictions with the stan- 
dard and the fine grid in the FPRBVI computation is 
presented. The large cylinder was used. The experi- 
mental data are also shown in the figures for reference. 
Although the fine-grid inputs yield slightly better re- 
sults (the difference is approximately 3 percent), the 
computing cost for the fine-grid solution with FPRBVI 
is substantially greater. Both codes predict, the magni- 
tude of the negative pressure peak with approximately 
the same accuracy. The agreement with the data is 
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Figure 2. Comparison of standard- and fine-grid 
FPRBVI inputs to WOPWOP+ (large cylinder) for 
UH-1H model rotor operation at Mh = 0.88. 
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Figure 3. Comparison of standard- and fine-grid 
FPRBVI inputs to RIvIR (large cylinder) for UH-1H 
model rotor operation at Mh = 0.88. 
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time 


Figure 4. Comparison of small- and large-cylinder in- 
puts to WOPWOP+ (fine grid) for UH-1H model rotor 
operation at Mh = 0.88. 

quite good for both methods. 

In Figs. 4 and 5, the acoustic predictions made by 
using the two surfaces shown in Fig. 1 are compared. 
The fine-grid results were utilized for this comparison 
because the fine grid yields slightly better agreement 
with the data. The tip-end surface was chosen to be 
just inside; the sonic cylinder for this Mach number. 
The results were relatively insensitive to the choice of 
the position of this end surface, as long as it was lo- 
cated near the sonic cylinder for high-speed cases. In 
Fig. 4, the WOPWOP+ calculation indicated nearly 
no change in predicted signal due to the increase in in- 
tegration volume. The reason that the predicted sig- 
nal does not change is apparently because the main 
qua.drupole contributions come primarily from the re- 
gion near the blade. The slight change in the negative 
peak pressure when the larger volume is used corre- 
sponds to the inclusion of more qua.drupole sources. 
The results for the Kirchhoff method in Fig. 5, how- 
ever, show a large sensitivity to the location of the 
Kirchhoff surface. These results show a nearly con- 
stant. shift, in the curves; this shift, is nearly identical 
for both the standard and the fine grid. One pos- 
sible explanation for this shift, is the nonlinearity of 
the input, signal. Another explanation may be that, 
the spatial pressure derivatives are not. computed with 
sufficient, accuracy on either grid. An important, goal 
of future work will be the determination of the cause 
of this variability. 
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Figure 5. Comparison of small- and large-cylinder in- 
puts to RKIR (fine grid) for UH-1H model rotor oper- 
ating at Mjj = 0.88. 

Now that we have examined the sensitivity of the 
codes for a representative case, we turn our atten- 
tion to comparison of method performance for a range 
of operating conditions. For the following cases, we 
present only predictions from the large cylinder on 
the standard CFD grid. These comparisons are repre- 
sentative of computations that can be routinely made 
with the FPRBVI solver. 

Hover Results 

For Mh = 0.70, the flow field around the rotor is 
subcritical; thus, the quadrupole contribution in the 
WOPWOP+ solution is not expected to be very large. 
Figure 6 shows the WOPWOP+ and RKIR predic- 
tions of the acoustic pressure •compared with the Euler 
calculation (TURNS code [3]) of Baeder et. al. [19]. 
The quadrupole component of the total signal pre- 
dicted with WOPWOP+ is also shown in Fig. 6. 
Notice that although the quadrupole contribution is 
small, the contribution is significant even for this low 
subcritical Mach number and should be included in an 
accurate acoustic prediction when the acoustic anal- 
ogy method is used. The Ivirchhoff prediction is nearly 
the same as the WOPWOP+ prediction. The excel- 
lent. agreement between the three predictions suggest 
that the FPRBVI solution is accurate. (The acous- 
tic analogy signal is dominated by thickness noise - 
which is computed from a knowledge of only the blade 
kinematics and geometry. Thus if the FPRBVI solu- 



time 

Figure 6. Comparison of acoustic predictions for UH- 
1H model rotor operation at Mh = 0.70 with the 
standard FPRBVI grid. (Reference computation from 
Baeder et al. [19].) 

t.ion is not sufficiently accurate at this Mach number, 
it will have only a small effect on the acoustic analogy 
prediction, but the Ivirchhoff solution will suffer.) 

The case for which Mh = 0.88 has already been 
examined, but we would like to show a direct com- 
parison of the WOPWOP+ and RKIR predictions. 
(See Fig. 7.) In this case, the quadrupole compo- 
nent. of the WOPWOP+ prediction is a. much higher 
percentage of the total (nearly 50 percent.) than for 
the cases in which Mh = 0.70. The Ivirchhoff predic- 
tion is nearly indistinguishable from the WOPWOP+ 
predict.ion. The slight, underpredict. ion of the nega- 
tive peak, as compared to the experimental data, is at. 
least, partially an indication that, the grid resolution of 
the standard grid is only marginally sufficient.. This 
inadequacy of the standard grid was also evident, in 
the earlier comparison of the standard- and fine-grid 
results. 

For M = 0.90, the flow becomes delocalized (i.e. , 
the supersonic region on the blade surface extends to 
a. region in which the free stream is moving super- 
sonically relative to the rotor-blade tip [8] ) . When 
delocalization occurs, shock waves on the rotor sur- 
face are allowed to propagate to the far field and the 
radiated noise dramatically increases. Because both 
methods are currently limited to subsonica.lly moving 
sources, the accuracy of t.hf yesult.s was expected t.o 
degrade in this case. Figure 8 shows the results for 
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Figure 7. Comparison of acoustic predictions for UH- 
1H model rotor operation at Mh = 0.88. 


Figure 8. Comparison of acoustic predictions for UH- 
1H model rotor operation at Mh = 0.90. 


Mh = 0.90 with the standard grid. For this case, 
the quadrupole component, is the largest component 
of the noise. Notice the dramatic increase in the peak 
negative pressure and the sawtooth waveform shape 
that results with a very small increase in Mach num- 
ber; these changes are characteristics of delocalization. 
The WOPWOP+ and RIvIR predictions for this con- 
dition are again in close agreement. Both codes fail 
to fully predict the sharp nonlinear steepening of the 
wave pulse. The negative peak the quadrupole contri- 
bution has the shape of the measured signal; however, 
the quadrupole component is not large enough. This 
underprediction results from not including the acous- 
tic sources that are moving supersonically. 

Advancing Rotor Results 

Two advancing rotor cases, chosen from separate ex- 
periments, are now described in this section. 

OLS Rotor 

The first, comparison for an advancing rotor is made 
for the Operational Loads Survey (OLS) model rotor. 
This rotor is a. one-seventh scale model of the main ro- 
tor of the AH-1 series helicopter. Model rotor acoustic 
data, were obtained by Schmitz et. a.l. [25] in the Duit.s- 
Nederla.ndse Windt.unnel (DNW) for several high for- 
ward speeds. Several microphones were used in the 
test., but. we will use data, from three in-pla.ne posi- 
tions located 1.72 rotor diameters from the rotor hub 
ahead of the rotor. The three microphones selected 
for use were in the following positions: 30° on the ad- 


vancing side of the rotor, on the centerline, and 30° on 
the retreating side of the rotor. The OLS rotor was 
a. t.wo-bla.ded rotor. For the following comparison, the 
rotor was operating with an advancing tip Mach num- 
ber Mat = 0.84 and an advance ratio p. = 0.27. A 
complete description of the test, can be found in Ref. 
[26]. 

Figures 9-11 show the results for this case at. the 
a.dva.ncing-side, centerline, and ret.rea.t.ing-side micro- 
phone positions, respectively. In all cases, we see that 
the results from both codes compare well with the ex- 
perimental measurements. The Ivirchhoff results are 
closer to experimental data, for all three microphone 
positions. Because peak negative pressure for this con- 
dition is not. large, the inclusion of loading in the pre- 
dictions may improve the waveform shape. 

UTC Four-bladed rotor 

The second advancing rotor case is for a. four-bla.ded 
United Technologies Corporation contemporary design 
model rotor. This rotor was pressure instrumented and 
tested in the DNW tunnel. The rotor had swept, tips 
and modern airfoil sections. The test, is described by 
Yu et. a.l. [27]; the acoustics portion of the test, is sum- 
marized by Liu and Ma.rcolini [28]. The measurements 
used in our comparison were taken in the tip-path 
plane of the rotor, at. 30° on the advancing side, on 
the centerline, and at. 30° on the retreating side ahead 
of the rotor, 1.5 diameters from the rotor hub. Com- 
parisons between the predicted and experimental data. 
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Figure 9. Comparison of advancing-side microphone 
predictions (standard grid) for OLS model rotor oper- 
ation at Mat = 0.84, and p, = 0.27. 




Figure 11. Comparison of retreating-side microphone 
predictions (standard grid) for OLS model rotor oper- 
ation at Mat = 0.84, and p, = 0.27. 


were made for a high-speed condition Mat = 0.94, 
and p. = 0.32. The flow around the rotor was delo- 
calized. A degraded accuracy was expected for this 
case because neither code accounts for the supersoni- 
cally moving sources; however, this case was expected 
to provide a good indication of code robustness. 

Figures 12-14 show the results for this case at each 
microphone location. The quality of the results over- 
all is good. Both codes predict lower values for the 
peak pressures than those measured experimentally, 
which can be attributed to the absence of the superson- 
ically moving sources. Nevertheless, the predictions 
are significantly better than would have been obtained 
by the WOPWOP code, which includes only thickness 
and loading sources. Note that the RIvIR prediction is 
better for the retreating-side microphone, whereas the 
WOPWOP+ prediction is better for the advancing- 
side microphone. The reason for this difference is un- 
clear. Thft Kirchhoff prediction does not agree well 
with either WOPWOP+ or the data away from the 
negative peak; this may indicate that the placement 
of the Kirchhoff surface is not optimal. 


Figure 10. Comparison of centerline microphone pre- 
dictions (standard grid) for OLS model rotor operation 
at Mat = 0.84, and p, = 0.27. 


Concluding Remarks 

Two methods were compared for computing ro- 
torcraft impulsive noise: an acoustic analogy ap- 

proach that includes the Ffowcs Williams-Hawkings 
quadrupole terms and a rotating Kirchhoff method. 
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Figure 12. Comparison of advancing-side microphone 
predictions (standard grid) for UTC model rotor op- 
eration at Mat = 0.93, and p, = 0.32. 
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Figure 13. Comparison of centerline microphone pre- 
dictions (standard grid) for UTC model rotor opera- 
tion at Mat = 0.93, and p, = 0.32. 
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Figure 14. Comparison of retreating-side microphone 
predictions (standard grid) for UTC model rotor op- 
eration at Mat = 0.93, and p, = 0.32. 


Both methods rely on computational fluid dynam- 
ics input. Because the CFD code and both acous- 
tic prediction codes are implemented in the Tilt.Ro- 
t.or Aeroacoustic Code (TRAC) at NASA Langley Re- 
search Center, the results in this paper should be use- 
ful to future users of this code system. Although the 
methods were coupled to a full-potential rotor finite- 
difference method, these methods can easily be cou- 
pled to other CFD solvers, as well. 

Both methods agree well with the experimental mea- 
surements for a variety of rotors and transonic operat- 
ing conditions below the delocalization Mach number. 
The Kirchhoff method is more sensitive to the choice 
of control surface; however, the choice of the Kirchhoff 
surface is nearly independent of the Mach number in 
the cases considered here. Because neither code in- 
cludes supersonically moving sources, performance is 
degraded for delocalized cases. For both codes, results 
with a fine grid are more accurate than the results 
with a standard grid, but the accuracy improvements 
are not dramatic. The computer time required by the 
flow solver, however, is substantially greater for the 
fine grid. 

The computational resources required by both com- 
puter codes is roughly equivalent. The Kirchhoff code 
requires less computer time for a single observer than 
WOPWOP+ in forward-flight cases, but the time to 
compute hover predictions is approximately the same. 
The computational requirements of the two acoustics 
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prediction methods are minimal, however, in compar- 
ison with the resources required by the CFD solver. 
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